Nonlinear control method for micro-grid inverter with anti-disturbance

ABSTRACT

Nonlinear control method for the micro-grid inverter with anti-disturbance. By generating reference currents that satisfy specific active and reactive power command under various working conditions, and introducing a nonlinear control method based on Lyapunov function to control the inverter, fast and accurate tracking of the generated reference signals is realized. The method realizes effective decoupling control of active power and reactive power. The system has high dynamic response and good robustness. Besides, the control structure of the method is simple and easy to implement, and the synchronous control link and the additional voltage and current regulator are omitted. The method realizes fast and accurate power exchange and stable power transmission between the inverter and the grid in the micro-grid under various working conditions, and provides a guarantee for improving the energy management efficiency within the micro-grid.

CROSS-REFERENCE TO RELATED APPLICATION

The subject application claims priority on Chinese patent application CN201810206558.3 filed on Mar. 13, 2018 in China. The contents and subject matter of the Chinese priority application is incorporated herein by reference.

FIELD OF TECHNOLOGY

The present invention relates to an electronic power converter, particularly, a nonlinear control method for a micro-grid inverter with anti-disturbance.

BACKGROUND ART

The micro-grid is one of the key technologies of the smart grid that can effectively reduce the impact on the main network and users due to the dispersed location, diverse forms, and characteristics of distributed power sources, and realize effective integration and efficient utilization of distributed power sources. As a result, the micro-grid has received widespread attention in recent years. The operation of the micro-grid is flexible, as it can operate in parallel to large grid or in island.

The micro-grid contains a large amount of renewable energy, and most of the renewable energy is incorporated into the micro-grid through inverters, therefore, it is especially important to control inverters with advanced technology. When the micro-grid is in grid-connected operation state, each inverter needs to accurately transmit or receive active and reactive power according to power commands issued by the central controller of the micro-grid, and thereby meets requirements of specific power exchange between the micro-grid and large grid and improves energy management efficiency within the micro-grid.

The micro-grid of medium- and low-voltage grades, however, is often at the end of the distribution network. The voltage quality of the common connection with the distribution network is not ideal. There are often numerous background harmonics, which seriously affects the power transmission quality of the grid-connected inverters. At the same time, sag and swell of the grid voltage amplitude, frequency fluctuations, and line impedance changes cause the inverter transmission power transmitted by the inverter to fluctuate, which deviates from the preset power command, resulting in the decrease in energy transfer efficiency and even leading to negative effect on the system stability and inverter integration.

At present, the control of the transmission power of the grid-connected inverter is mainly realized by controlling the output current of the inverter. The traditional control method is mainly double-loop control. The outer loop controller generates the reference signal. The inner loop is mostly the current loop, which is related to the stability precision, harmonic content, dynamic response, and anti-disturbance ability of the inverter. According to different coordinate system selection, different control strategies are needed. For example, in the synchronous rotating coordinate system, AC variables need to be converted into DC variables, and the proportional integral controller (PI) is used to eliminate the steady-state error. However, in the non-ideal state of the grid voltage, variables obtained by the synchronous rotating coordinate system conversion is no longer DC variables, and traditional PI control is not able to meet requirements of zero steady-state error tracking. And the control in the synchronous rotating coordinate system involves phase-locking part. The performance of phase-locking part directly affects the dynamic response speed of the system and the control effect of the inverter. If the control is performed in the two-phase static coordinate system, the synchronization part required for the transformation can be eliminated, and the complexity of the control system can be reduced. The proportional integral regulator (PR) can exhibit high gain characteristics at the fundamental frequency, which can approximatively achieve error tracking for a sinusoidal quantity of a specific frequency. However, when the sinusoidal component contains other harmonic components, the resonance controller has to be incorporated at specific frequency to eliminate other harmonic components, which increases the complexity of the controller. In the control under the three-phase natural coordinate system, the three-phase system is divided into three single-phase systems for control, thus the control of each phase is relatively independent. Although it can be used for three-phase unbalanced systems, the structure of the controller is relatively complicated.

At the same time, with the in-depth study of nonlinear control theory, some nonlinear control performance characteristics such as fast dynamic response, global stability, and strong robustness, etc., to some extent, make up for the lack of linear control. In recent years, the characteristics have attracted more and more scholars' attention, and have also been applied in the control of inverters, but there are still many problems need to be solved.

SUMMARY OF THE PRESENT INVENTION

In order to overcome the deficiency of the traditional methods, the present invention provides a nonlinear control method for the micro-grid with anti-disturbance to improve the ability of the grid-connected inverter to cope with the poor power quality environment at the grid-connected point, so that each inverter is able to emit active and reactive power steadily and accurately according to the issued power command in the presence of interference conditions, thereby improving the reliability of the inverter operation and the operating efficiency of the micro-grid. The present invention provides a reasonably designed nonlinear control method for power of the LCL type grid-connected inverter with good dynamic and steady state characteristics.

To achieve the objective of the present invention, the present invention provides a nonlinear control method for the micro-grid inverter with anti-disturbance, which includes following steps:

Step 1: Collect the LCL type grid-connected inverter filter capacitor voltage v_(c) and the grid-side inductance current i_(o), and the mathematical model of the two parameters in the two-phase stationary coordinate system are established through coordinate transformation.

Step 2: Collect the AC bus voltage v_(g) and the local load current i_(load). After coordinate transformation, the grid-side inductance current reference signal in the two-phase stationary coordinate system is generated according to the instantaneous reactive power theory. The derivative of the reference signal is taken by the observer to obtain the first and second derivative value of the reference signal as input signal values of the nonlinear controller.

Step 3: Set an original tracking error and a filter tracking error of the control variable i_(o) to establish a nonlinear control model of the grid-connected inverter.

Step 4: The modulated wave signal is obtained according to the LCL type grid-connected inverter mathematical model and the grid-connected inverter nonlinear control model. After coordinate transformation, the modulated wave signal in the three-phase stationary coordinate system is obtained. The grid-connected inverter is controlled by introducing a nonlinear control method based on Lyapunov function.

The mathematical model of the grid-connected inverter in a two-phase stationary coordinate system is as follows:

$\begin{matrix} {{\frac{{di}_{o\; \alpha \; \beta}}{dt} = {x_{\alpha \; \beta} - d}};{\frac{{dx}_{\alpha \; \beta}}{dt} = u_{\alpha \; \beta}}} & (1) \end{matrix}$

The method for generating the current reference signal in the two-phase stationary coordinate system according to the instantaneous reactive power theory is as follows:

$\begin{matrix} {i_{o\; \alpha \; \beta \; \_ \; {ref}} = {\begin{bmatrix} i_{o\; \alpha \; \_ \; {ref}} \\ i_{o\; \beta \; \_ \; {ref}} \end{bmatrix} = {\begin{bmatrix} i_{{load}\mspace{14mu} \alpha} \\ i_{{load}\mspace{14mu} \beta} \end{bmatrix} - {\frac{2}{3} \cdot {{\frac{1}{v_{g\; \alpha}^{2} + v_{v\; \beta}^{2}}\begin{bmatrix} v_{g\; \alpha} & v_{g\; \beta} \\ v_{g\; \beta} & {- v_{g\; \alpha}} \end{bmatrix}}\begin{bmatrix} P_{set} \\ Q_{set} \end{bmatrix}}}}}} & (2) \end{matrix}$

The original tracking error and the filter tracking error are respectively set as follows:

eαβ=i _(oαβ) −i _(oαβ_ref) ;E _(αβ) =ke _(αβ) +ė _(αβ)  (3)

The grid-connected inverter nonlinear control model is as follows:

$\begin{matrix} {{u_{\alpha \; \beta}(t)} = {{- {{kx}_{\alpha \; \beta}(t)}} - {\left( {\lambda + k_{s}} \right){E_{\alpha \; \beta}(t)}} + {k{{\overset{.}{i}}_{o\; \alpha \; \beta \; \_ \; {ref}}(t)}} + {{\overset{¨}{i}}_{o\; \alpha \; \beta \; \_ \; {ref}}(t)} - {\int_{0}^{t}{\left( {{\lambda \; k_{s}} + \frac{1}{\rho}} \right){E_{\alpha \; \beta}(\tau)}d\; \tau}}}} & (4) \end{matrix}$

where i_(oαβ) represents the grid-side inductance current component in the two-phase stationary coordinate system; x_(αβ) represents the state variable, d is the disturbance, including the error caused by the system disturbance and the model uncertainty; i_(loadαβ) represents the load current component in the two-phase stationary coordinate system; v_(gαβ) represents the grid voltage component in the two-phase stationary coordinate system; P_(set), Q_(set) are the given active power and reactive power commands respectively; i_(oαβ_ref) is the reference current signal component; e_(αβ) and E_(αβ) are the set original error and filter tracking error respectively; u_(αβ) represents the nonlinear control law in the two-phase stationary coordinate system, where k, λ, k_(s), p are custom constants greater than zero, respectively.

The state variable is expressed as follows:

x _(αβ) =v _(cαβ) /L _(o)  (5)

The disturbance is as follows:

$\begin{matrix} {d = {{\frac{1}{L_{o}}\left( {v_{g\; \alpha \; \beta} + {\Delta \; v_{g\; {\alpha\beta}}}} \right)} + {\frac{i_{o\; {\alpha\beta}}}{L_{o}}\left( {R_{o} + {\Delta \; R_{o}}} \right)} + {\frac{1}{L_{o}}\left( {{\Delta \; L_{o}{\overset{.}{i}}_{o\; \alpha \; \beta}} - {\Delta \; v_{c\; \alpha \; \beta}}} \right)}}} & (6) \end{matrix}$

where L _(o) is an estimate of the sum L_(o) of the grid-side filter inductance and the line inductance; Δv_(gαβ) is the disturbance caused by the grid voltage; ΔR_(o) and ΔL_(o) are the deviations between the actual applied value and the theoretical value due to the uncertainty of the line parameters; Δv_(cαβ) is the effect of the grid voltage disturbance on the capacitor voltage.

The generated modulated wave signal is as follows:

$\begin{matrix} {v_{s\; \alpha \; \beta} = {{\frac{1}{K_{d\;}} \cdot v_{i\; \alpha \; \beta}} = {\frac{1}{K_{d}} \cdot \left\{ {{\frac{d}{dt}\left( {{L_{o}{C \cdot u_{\alpha \; \beta}}} + i_{o\; \alpha \; \beta}} \right)L_{i}} + v_{c\; \alpha \; \beta}} \right\}}}} & (7) \end{matrix}$

where v_(sαβ) represents the modulated wave signal component in the two-phase stationary coordinate system; v_(iαβ) represents the inverter input voltage signal component in the two-phase stationary coordinate system; K_(d) is the scaling factor and is greater than zero; L _(i), C, L _(o) are estimates of the inverter-side filter inductance L_(i), the filter capacitor C, and the sum L_(i) of grid-side filter inductance and the line parameter; v_(cαβ) represents the capacitor voltage component in the two-phase stationary coordinate system.

By reasonable value of the custom parameters in the nonlinear control model, the stability of the system is proved by the Lyapunov's direct method and the tracking error tends to zero.

Lyapunov function: V=½E²+½ρr² is selected and the value of its derivative {dot over (V)} is related to custom parameters k, λ, k_(s) and ρ. By the selection of custom parameters, the derivative {dot over (V)} of V is less than zero. As V is a function of E and r and greater than zero, its derivative {dot over (V)} is less than zero to ensure that V monotonically decreases until E and r tend to zero. Therefore, the global stability of the system is proved, and the tracking error tending to be zero is realized.

The disturbance d includes modeling of the uncertainty of the line impedance parameters ΔR_(o) and ΔL_(o), and ensures the power control accuracy of the designed controller when the system parameters are uncertain or the line impedance changes; the modeling of the disturbance on grid Δv_(g) and the influence of grid disturbance on the capacitor voltage Δv_(c) is included to ensure the control performance of the designed controller under non-ideal grid voltage conditions.

The present invention is advantageous in that:

1. It considers the interference on the LCL type grid-connected inverter caused by various factors such as poor voltage quality (more harmonics), voltage amplitude fluctuation, and frequency hopping at the grid-connected point, and the stability of the system is proved by the Lyapunov's direct method. In the present invention, the inverter is ensured to transmit or receive active and reactive power, the preset value is accurately tracked, the power fluctuation caused by the grid voltage unbalance is suppressed, and the anti-interference ability of the inverter is improved.

2. The present invention is realized in a two-phase stationary coordinate system, and the phase locking part is omitted, so that the system has better dynamic performance. When the active or reactive power command saltation occurs, the system can respond quickly, so as to achieve fast tracking of power, and ensure the transient stability of the system.

3. When the load is nonlinear, the quality of the grid-connected current is ensured, and the accurate transmission of active or reactive power is realized. Equations (4)-(5) are formulas for reference voltage. When the load is nonlinear, the load current waveform is non-sinusoidal. To ensure that the grid-connected current is a standard sine wave, the obtained reference value of the inverter output current must be non-sinusoidal. That is, i_(oαβ_ref)=i_(loadαβ)−i_(gαβ). The designed controller accurately tracks arbitrary waveforms, thus ensuring the quality of the grid-connected current.

4. The present invention considers parameter uncertainties of the LCL type grid-connected inverter system. When the line impedance changes, the system has strong robustness and the accuracy of power transmission are not affected.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic diagram of a micro-grid system including an LCL type grid-connected inverter of the present invention.

FIG. 2 is a schematic diagram showing the control method of the present invention.

FIGS. 3(a) to 3(c) show the voltage waveform and the system active/reactive waveforms when the grid-connected point is interfered according to the present invention, among which, FIG. 3(a) shows the voltage waveform and the system active/reactive waveforms when harmonic is contained in the voltage of the grid-connected point; FIG. 3(b) shows the voltage waveform and the system active/reactive waveforms when voltage amplitude of the grid connection increases/decreases; and FIG. 3(c) shows the voltage waveform and the system active/reactive waveforms when frequency hopping of the grid-connected point occurs.

FIG. 4 shows waveforms of the system active/reactive power when the active/reactive power command changes suddenly in the present invention, indicating that the active/reactive power of the system has fast dynamic responses without any delay and transients to keep up with changes in power command.

FIG. 5 shows the active/reactive power and current waveforms of the system under nonlinear loads in the present invention, indicating that the active/reactive power of the system has little affected by load changes even the loads are nonlinear.

DETAILED DESCRIPTION OF THE PRESENT INVENTION AND EMBODIMENTS

In combination with figures, the embodiments of the present invention is further expounded.

The nonlinear control method for the micro-grid inverter with anti-disturbance of the present invention is implemented on an LCL type grid-connected inverter in the micro-grid system shown in FIG. 1. The system mainly includes six parts: an external power grid 1, a grid-connected point switch 2, a three-phase AC bus 3, a three-phase load 4, a distributed electric generation element 5, and other distributed electric generation elements 6.

The distributed electric generation element 5 includes a DC voltage source 7, a three-phase voltage source full-bridge inverter 8, an LCL-type filter circuit 9, and a line impedance 10. The DC voltage source 7 uses an ideal DC voltage source, and the voltage is expressed as V_(dc); the three-phase voltage source full-bridge inverter 8 includes six switch tubes S₁˜S₆; the LCL-type filter circuit 9 includes an inverter-side filter inductance L_(i) 11, a filter capacitor C 12, and a grid-side filter inductance L_(g) 13; the line impedance 10 includes the line resistance R_(o) 14 and the line inductance L_(l) 15. For the modeling's sake, the grid-side filter inductance L_(g) and the equivalent line inductance L_(l) are combined into L_(o) 16. The no-load voltage of the three-phase voltage source full-bridge inverter 8 is v_(i); the voltage of the filter capacitor 12 is v_(c); the voltage at the grid-connected point is v_(g), the inductance current on the inverter side is i_(i); the inductance current on the grid side is i_(o); the load current is i_(load); the current flowing into the grid is i_(g).

The control method of the present invention includes the following steps:

Step 1: Collect the LCL type grid-connected inverter filter capacitor voltage v_(c) and the grid-side inductance current i_(o), and the mathematical model of the two parameters in the two-phase stationary coordinate system are established through coordinate transformation;

First, according to FIGS. 1 and 2, the mathematical model of the system in the two-phase stationary coordinate system is expressed as follows:

$\begin{matrix} {{{L_{i}\frac{{di}_{i\; {\alpha\beta}}}{dt}} = {v_{i\; {\alpha\beta}} - v_{c\; \alpha \; \beta}}}{{L_{o}\frac{{di}_{o\; \alpha \; \beta}}{dt}} = {v_{c\; \alpha \; \beta} - v_{g\; \alpha \; \beta} - {R_{o}i_{o\; \alpha \; \beta}}}}{{C\; \frac{{dv}_{c\; \alpha \; \beta}}{dt}} = {i_{i\; \alpha \; \beta} - i_{o\; {\alpha\beta}}}}} & (8) \end{matrix}$

Considering the uncertainty of model parameters and the interference that the system may experience, equation (17) can be expressed as follows:

$\begin{matrix} {{\overset{.}{i}}_{o\; {\alpha\beta}} = {\frac{v_{c\; {\alpha\beta}}}{L_{o}} - {\frac{1}{L_{o}}\left( {v_{g\; {\alpha\beta}} + {R_{o}i_{o\; {\alpha\beta}}}} \right)} - {\frac{1}{L_{o}}\left( {{\Delta \; L_{o}{\overset{.}{i}}_{o\; {\alpha\beta}}} + {\Delta \; v_{g\; {\alpha\beta}}} + {\Delta \; R_{o}i_{o\; {\alpha\beta}}} - {\Delta \; v_{c\; {\alpha\beta}}}} \right)}}} & (18) \end{matrix}$

where ΔL_(o), ΔR_(o) represents the deviation of system parameter estimation; Δv_(gαβ) represents the disturbance caused by the grid voltage, and Δv_(cαβ) represents the capacitor voltage disturbance caused by system parameter and grid voltage.

The disturbance is expressed as follows:

$\begin{matrix} {{{d_{1} = {\frac{1}{L_{o}}\left( {v_{g\; {\alpha\beta}} + {\Delta \; v_{g\; {\alpha\beta}}}} \right)}};{d_{2} = {\frac{i_{o\; {\alpha\beta}}}{L_{o}}\left( {R_{o} + {\Delta \; R_{o}}} \right)}};}d_{3} = {\frac{1}{L_{o}}\left( {{\Delta \; L_{o}{\overset{.}{i}}_{o\; {\alpha\beta}}} - {\Delta \; v_{c\; {\alpha\beta}}}} \right)}} & (9) \end{matrix}$

Define state variables:

${x_{\alpha\beta} = \frac{v_{c\; {\alpha\beta}}}{{\overset{\_}{L}}_{o}}},$

the control law is expressed as:

${u_{\alpha\beta} = {\frac{1}{{\overset{\_}{L}}_{o}\overset{\_}{C}}\left( {i_{i\; {\alpha\beta}} - i_{o\; {\alpha\beta}}} \right)}},$

where L _(o), C are estimates of L_(o), C, respectively. Therefore, the simplified mathematical model of the grid-connected inverter is as follows:

$\begin{matrix} {{\frac{{di}_{o\; {\alpha\beta}}}{dt} = {x_{\alpha\beta} - d}};{\frac{{dx}_{\alpha\beta}}{dt} = u_{\alpha\beta}}} & (10) \end{matrix}$

where d is the disturbance and can be expressed as

${d = \frac{d_{1} + d_{2} + d_{3}}{3}},$

which is assumed to be unknown and finite.

Step 2: Collect the AC bus voltage v_(g) and the local load current i_(load). After coordinate transformation, the grid-side inductance current reference signal in the two-phase stationary coordinate system is generated according to the instantaneous reactive power theory. The derivative of the reference signal is taken by the observer to obtain the first and second derivative value of the reference signal as input signal values of the nonlinear controller;

The collected signals v_(g) and i_(load) are coordinate transformed to obtain v_(gαβ) and i_(loadαβ), and the current reference value i_(oαβ_ref) of the grid-side filter inductance is obtained by the calculation formula (21) based on the instantaneous reactive power theory.

$\begin{matrix} {i_{o\; {{\alpha\beta\_}{ref}}} = {\begin{bmatrix} i_{o\; {{\alpha\_}{ref}}} \\ i_{o\; {{\beta\_}{ref}}} \end{bmatrix} = {\begin{bmatrix} i_{{load}\; \alpha} \\ i_{{load}\; \beta} \end{bmatrix} - {\frac{2}{3} \cdot {{\frac{1}{v_{g\; \alpha}^{2} + v_{g\; \beta}^{2}}\begin{bmatrix} v_{g\; \alpha} & v_{g\; \beta} \\ v_{g\; \beta} & {- v_{g\; \alpha}} \end{bmatrix}}\begin{bmatrix} P_{set} \\ Q_{set} \end{bmatrix}}}}}} & (11) \end{matrix}$

Then, accurate estimates of the first and second derivatives of the reference signal are obtained by the observer:

{circumflex over ({dot over (i)})}_(oαβ_ref) =k _(o)(i _(oαβ_ref) −î _(oαβ_ref))  (12)

where {circumflex over ({dot over (i)})}_(oαβ_ref) and î_(oαβ_ref) are the estimated values of {dot over (i)}_(oαβ_ref) and i_(oαβ_ref), respectively, k_(o) is the gain of the observer, and î_(oαβ_ref) is updated in real time according to {dot over (î)}_(oαβ_ref)=k_(o)(i_(oαβ_ref)−î_(oαβ_ref)).

Define the output error of the observer as follows:

e ₁ =i _(oαβ_ref) −î _(oαβ_ref) ;e ₂ ={dot over (i)} _(oαβ_ref)−{circumflex over ({dot over (i)})}_(oαβ_ref)  (13)

By setting the gain of the observer k₀ to a large enough value, the output error of the observer can be adjusted to within the range that the system allows. In the same way, the second derivative of the current reference signal is obtained on the basis of the first derivative.

Step 3: Set the original tracking error and the filter tracking error of the control variable i_(o) to establish a nonlinear control model of the grid-connected inverter;

The original error signal and the filter tracking error signal are defined as:

e _(αβ) =i _(oαβ) −i _(oαβ_ref) ;E _(αβ) =ke _(αβ) +ė _(αβ)  (14)

Through the observer in step 2, the first derivative {dot over (i)}_(oαβ) of the grid-side filter inductance current is obtained. The filter tracking error E_(αβ) is obtained as an input signal of the nonlinear controller by the formula (24), and the nonlinear control law is as follows:

$\begin{matrix} {{u_{\alpha\beta}(t)} = {{- {{kx}_{\alpha\beta}(t)}} - {\left( {\lambda + k_{s}} \right){E_{\alpha\beta}(t)}} + {k{{\overset{.}{i}}_{o\; {{\alpha\beta\_}{ref}}}(t)}} + {{\overset{¨}{i}}_{o\; {{\alpha\beta\_}{ref}}}(t)} - {\int_{0}^{t}{\left( {{\lambda \; k_{s}} + \frac{1}{\rho}} \right){E_{\alpha\beta}(\tau)}d\; \tau}}}} & (15) \end{matrix}$

where k, λ, k_(s), and ρ are custom constants greater than zero, respectively.

Step 4: The modulated wave signal is obtained according to the LCL type grid-connected inverter mathematical model and the grid-connected inverter nonlinear control model. After coordinate transformation, the modulated wave signal in the three-phase stationary coordinate system is obtained. The grid-connected inverter is controlled by introducing a nonlinear control method based on Lyapunov function.

According to equation (17), the modulation wave signal in the two-phase stationary coordinate system is as follows:

$\begin{matrix} {v_{s\; {\alpha\beta}} = {{\frac{1}{K_{d}} \cdot v_{i\; {\alpha\beta}}} = {\frac{1}{K_{d}} \cdot \left\{ {{\frac{d}{dt}\left( {{{\overset{\_}{L}}_{o}{\overset{\_}{C} \cdot u_{\alpha\beta}}} + i_{o\; {\alpha\beta}}} \right){\overset{\_}{L}}_{i}} + v_{c\; {\alpha\beta}}} \right\}}}} & (16) \end{matrix}$

where L _(i) is the estimated value of the inverter-side filter inductance; K_(d) is the scaling factor greater than zero.

The modulation wave signal is input into the sine wave vector modulation module to obtain pulse signals of six channels, S₁˜S₆, which are input into the three-phase voltage source full-bridge inverter circuit to control the said circuit.

In conclusion, the nonlinear power control method of the LCL type grid-connected inverter can be obtained, as shown in FIG. 2.

A simulation model of the LCL type grid-connected inverter system is built by Matlab/Simulink.

Three working conditions are simulated to test the control method of the present invention, that is, the simulated inverter is affected by the voltage disturbance (containing harmonics, voltage amplitude fluctuation, and frequency hopping) at the grid-connected point, the power preset value is changed, and the nonlinear load is connected.

Working condition 1: The grid voltage at the grid-connected point is disturbed. There are three main situations:

1. Background harmonics are contained;

2. Voltage amplitude fluctuation;

3. Frequency fluctuation.

The power value is set to input 20 kW of active power to the grid and transmit 10 kVar of reactive power. The effects of the above three grid-side disturbances on the power input to the grid by the inverter are verified. The simulation results are shown in FIGS. 3(a) to 3(c). v_(ga), v_(gb), and v_(gc) represent voltages of three phases at the grid-connected point abc respectively. P_(g) and Q_(g) represent the active and reactive power exchange between the inverter and the grid respectively. The simulation time of the system is 0.2 s. For the convenience of comparing, the voltage at the grid-connected point during 0˜0.1 s is a ideal voltage waveform, and the voltage at the grid-connected point during 0.1 s˜0.2 s is disturbed, which is a non-ideal voltage waveform.

Firstly, after 0.1 s, it is set that the voltage at the grid-connected point contains the third harmonic component of 15V and the fifth harmonic component of 15V, while the phase of the third harmonic component phase lags the phase of the fundamental component by 25 degrees, and the phase of the fifth harmonic component phase leads the phase of the fundamental component by 35 degrees. It can be seen from FIG. 3(a) that after 0.1 s, the active and reactive power input into the grid by the inverter can accurately track the preset power command, and the power transmission is kept stable when the voltage at the grid-connected point is distorted. No second harmonic pulsation is generated and no large fluctuation occurs at the moment when the voltage condition deteriorates.

Then, after 0.1 s, it is set that the voltage amplitude at the grid-connected point begin to fluctuate. It can also be observed that the active and reactive power input into the grid are always maintained around the preset value and remains stable at all times.

Finally, after 0.1 s, it is set that the voltage frequency at the grid-connected point jumps from 50 Hz to 49.5 Hz. It can be seen that the power transmission is not affected by the frequency change.

Working condition 2: The preset power command changes over time. In 0˜0.1 s, the active and reactive power input into the grid are set to zero; in 0.1 s˜0.2 s, the system transmits 10 kW of active power to the grid and receives 20 kVar of reactive power from the grid; in 0.2 s˜0.3 s, the system transmits 20 kW of active power to the grid and receives 10 kVar of reactive power from the grid; in 0.3 s˜0.4 s, the system transmits 10 kW of active power to the grid and receives 10 kVar of reactive power from the grid.

It can be seen from the simulation results that the system realizes decoupling control of active and reactive power, and can track power commands quickly and accurately.

Working condition 3: The system supplies power to the linear load from 0 to 0.1 s, and supplies power to the nonlinear load from 0.1 s to 0.2 s. And the power command is to transmit 20 kW of active power and 10 kVar of reactive power to the grid.

From the waveforms of the system active/reactive power when the active/reactive power command changes suddenly in the present invention as shown in FIG. 4, the active/reactive power of the system has fast dynamic responses without any delay and transients to keep up with changes in power command. From the active/reactive power and current waveforms of the system under nonlinear loads in the present invention as shown in FIG. 5, the active/reactive power of the system has little affected by load changes even the loads are nonlinear.

It can be seen from the simulation results that at the moment when the load sudden changes, the power has received a small impact and quickly converges to the command value, and ensures the smooth and accuracy quality of power transmission to the grid when the load is nonlinear.

It can be seen that by using the said control method, the LCL type grid-connected inverter has good power transmission capability and the dynamic response of the system is rapid under conditions that the voltage point at the grid-connected is disturbed or the load is nonlinear. 

We claim:
 1. A nonlinear control method for a micro-grid inverter with anti-disturbance, comprising: (1) collecting a filter capacitor voltage v_(c) of an LCL type grid-connected inverter and a grid-side inductance current i_(o), and building a mathematical model of the grid-connected inverter in a two-phase stationary coordinate system through coordinate transformation as follows: $\begin{matrix} {{\frac{{di}_{o\; {\alpha\beta}}}{dt} = {x_{\alpha\beta} - d}};{\frac{{dx}_{\alpha\beta}}{dt} = u_{\alpha\beta}}} & (17) \end{matrix}$ wherein i_(oαβ) represents a grid-side inductance current in the two-phase stationary coordinate system; u_(αβ) represents a nonlinear control law in the two-phase stationary coordinate system; x_(αβ) is a state variable as follows: x _(αβ) =v _(cαβ) /L _(o)  (18) d is a disturbance including error caused by system disturbance and model uncertainty as follows: $\begin{matrix} {d = {{\frac{1}{L_{o}}\left( {v_{g\; {\alpha\beta}} + {\Delta \; v_{g\; {\alpha\beta}}}} \right)} + {\frac{i_{o\; {\alpha\beta}}}{L_{o}}\left( {R_{o} + {\Delta \; R_{o}}} \right)} + {\frac{1}{L_{o}}\left( {{\Delta \; L_{o}{\overset{.}{i}}_{o\; {\alpha\beta}}} - {\Delta \; v_{c\; {\alpha\beta}}}} \right)}}} & (19) \end{matrix}$ wherein L _(o) is an estimate of a sum L_(o) of a grid-side filter inductance and a line inductance; Δv_(gαβ) is a disturbance caused by a grid voltage; R_(o) is a value of a line resistance; ΔR_(o) and ΔL_(o) are deviations between an actual applied value and a theoretical value due to uncertainty of line parameters; Δv_(cαβ) is an effect of a grid voltage disturbance on a capacitor voltage; i_(loadαβ) represents the load current in the two-phase stationary coordinate system; v_(gαβ) represents a grid voltage in the two-phase stationary coordinate system; (2) collecting an AC bus voltage v_(g) and a local load current i_(load) and generating a grid-side inductance current reference signal i_(oαβ_ref) in the two-phase stationary coordinate system according to an instantaneous reactive power theory after coordinate transformation as follows: $\begin{matrix} {{i_{o\; {{\alpha\beta\_}{ref}}} = {\begin{bmatrix} i_{o\; {{\alpha\_}{ref}}} \\ i_{o\; {{\beta\_}{ref}}} \end{bmatrix} = {\begin{bmatrix} i_{{load}\; \alpha} \\ i_{{load}\; \beta} \end{bmatrix} - \begin{bmatrix} i_{g\; \alpha} \\ i_{g\; \beta} \end{bmatrix}}}},} & (20) \end{matrix}$ calculating a grid-connected current i_(gαβ) by equation (5): $\begin{matrix} {i_{g\; {\alpha\beta}} = {\begin{bmatrix} i_{g\; \alpha} \\ i_{g\; \beta} \end{bmatrix} = {\frac{2}{3} \cdot {{\frac{1}{v_{g\; \alpha}^{2} + v_{g\; \beta}^{2}}\begin{bmatrix} v_{g\; \alpha} & v_{g\; \beta} \\ v_{g\; \beta} & {- v_{g\; \alpha}} \end{bmatrix}}\begin{bmatrix} P_{set} \\ Q_{set} \end{bmatrix}}}}} & (21) \end{matrix}$ wherein P_(set) and Q_(set) are given active power and reactive power commands, respectively; (3) setting an original tracking error e_(αβ) and a filter tracking error E_(αβ) of a control variable i_(o) as follows: e _(αβ) =i _(oαβ) −i _(oαβ_ref) ;E=E _(αβ) =ke _(αβ) +ė _(αβ)  (22) wherein i_(oαβ_ref) is the reference current signal; k is a custom constant greater than zero; ė_(αβ) is a first derivative of the original tracking error e_(αβ); defining a tracking error as follows: r=r _(αβ) =αE _(αβ) +Ė _(αβ)  (23) wherein a derivative of the reference current signal i_(oαβ_ref)(t) is taken by an observer to obtain a first derivative {dot over (i)}_(oαβ_ref)(t) and a second derivative ï_(oαβ_ref)(t) of the reference signal as input signal values of the nonlinear controller and build a grid-connected inverter nonlinear control model as follows: $\begin{matrix} {{u_{\alpha\beta}(t)} = {{- {{kx}_{\alpha\beta}(t)}} - {\left( {\lambda + k_{s}} \right){E_{\alpha\beta}(t)}} + {k{{\overset{.}{i}}_{o\; {{\alpha\beta\_}{ref}}}(t)}} + {{\overset{¨}{i}}_{o\; {{\alpha\beta\_}{ref}}}(t)} - {\int_{0}^{t}{\left( {{\lambda \; k_{s}} + \frac{1}{\rho}} \right){E_{\alpha\beta}(\tau)}d\; \tau}}}} & (24) \end{matrix}$ wherein i_(loadαβ) represents the load current in the two-phase stationary coordinate system; k, λ, k_(s), p are custom constants greater than zero, respectively; (4) obtaining a modulated wave signal according to the LCL type grid-connected inverter mathematical model and the grid-connected inverter nonlinear control model in a three-phase stationary coordinate system after coordinate transformation as follows: $\begin{matrix} {v_{s\; {\alpha\beta}} = {{\frac{1}{K_{d}} \cdot v_{i\; {\alpha\beta}}} = {\frac{1}{K_{d}} \cdot \left\{ {{\frac{d}{dt}\left( {{{\overset{\_}{L}}_{o}{\overset{\_}{C} \cdot u_{\alpha\beta}}} + i_{o\; {\alpha\beta}}} \right){\overset{\_}{L}}_{i}} + v_{c\; {\alpha\beta}}} \right\}}}} & (25) \end{matrix}$ wherein v_(iαβ) represents an inverter input voltage signal in the two-phase stationary coordinate system; K_(d) is a scaling factor that is greater than zero and a DC voltage value of an inverter side; L _(i), C, L _(o) are estimates of the inverter-side filter inductance L_(i), the filter capacitor C, and the sum L_(i) of grid-side filter inductance and the line parameter; v_(cαβ) represents the capacitor voltage component in the two-phase stationary coordinate system; and (5) controlling the grid-connected inverter by introducing a nonlinear control method based on Lyapunov function and using the modulation wave signal in the three-phase stationary coordinate system. 